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ABSTRACT 

According to the cosmological principle, Universal large-scale structure is homoge¬ 
neous and isotropic. The observable Universe, however, shows complex structures even 
on very large scales. The recent discoveries of structures significantly exceeding the 
transition scale of 370 Mpc pose a challenge to the cosmological principle. 

We report here the discovery of the largest regular formation in the observable 
Universe; a ring with a diameter of 1720 Mpc, displayed by 9 gamma ray bursts 
(GRBs), exceeding by a factor of five the transition scale to the homogeneous and 
isotropic distribution. The ring has a major diameter of 43° and a minor diameter of 
30° at a distance of 2770 Mpc in the 0.78 < z < 0.86 redshift range, with a probability 
of 2 X 10“® of being the result of a random fluctuation in the GRB count rate. 

Evidence suggests that this feature is the projection of a shell onto the plane of the 
sky. Voids and string-like formations are common outcomes of large-scale structure. 
However, these structures have maximum sizes of 150 Mpc, which are an order of 
magnitude smaller than the observed GRB ring diameter. Evidence in support of 
the shell interpretation requires that temporal information of the transient GRBs be 
included in the analysis. 

This ring-shaped feature is large enough to contradict the cosmological principle. 
The physical mechanism responsible for causing it is unknown. 

Key words: Large-scale structure of Universe, cosmology: observations, gamma-ray 
burst: general 


1 INTRODUCTION 


Quasars are well-suited for mapping out the large-scale 
distribution of matter in the Universe, due to their very 
high luminosities and preferentially large redshifts. Quasars 
are associated by groups and poor clus ters of galaxies 


llHeinamaki et ah l2005l : iLietzen et al.ll2()()9l l and can be ob¬ 
served even when the underlying galaxies are faint and dif¬ 
ficult to detect. When quasars cluster, they identify con¬ 
siderable amounts of underlying matter, such that quasar 
clusters have been used to detect matter clustered on very 
large scales. Some of this matter is clustered on scales equal 
to or exceeding that of the Sloan Great Wall jGott et al.l 

I2OO5II . 

A number of large quasar groups (LQG) have been iden- 
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tified in recent years; each o ne mapping out large amounts 
of much fainter matter. After IWebsteJ (Il982l l found a group 
of four quasars at z = 0.37 with a size of about 100 
Mpc, having a low probability of bei ng a chance alignment, 
iKomberg. Kravtsov fc LukashI l|l994ll identified strong clus¬ 
tering in the quasar distribution at scales less than 20 h~^ 
Mpc, and defined LQGs usi ng a well-known cluster analy¬ 
sis tec hnique. Subsequently, iKomberg. Krayteoy_jMliukagh 
(1996) identified additional LQGs. and IKomberg fc Lukash 
1 199811 reported a new finding of eleven LQGs based on sys¬ 
tematic cluster analysis. T he sizes of these cluster s ranged 
from 70 to 160 h~^ Mpc. I Newman et al.l (Il998al f9l later 
discovered a 150 h~^ Mpc group of 1 3 quasars at median 
redshift z = 1.51. IWilliger et al.l (l2002l i mapped 18 quasars 
spanning « 5° x 2.5° on the sky, with a qu asar spatial over¬ 
densit y 6 — 10 times greater than the mean. lHaberzettl et al.l 
koogli investigated two sheet-like structures of galaxies at 
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z — 0.8 and 1.3 spanning 150comoving Mpc em- 
bedded in LQGs extending o ver a t least 200 h~^ Mpc. 
iHaines. Campusano fc Clowes! (l2004l reported the finding 
of two large-scale structures of galaxies in a 40 x 35 arcmin^ 
field embedded in a 25°^ area containing two 100 Mpc- 
scale structures of quas ars. As the identified scales of quasar 
clusters became larger, IClowes et al.l ll2012l l found two rel¬ 
atively close LQGs at 2 : ~ 1.2. The characteristic sizes of 
these two LQGs, « 350 — 400 Mpc, and appear to be only 
marginally consistent with the sc ale of homogeneity in the 
concordance cosmology. Recently, iGlowes et al.l (l2013l l un¬ 
covered the Huge-LQG with a long dimension of « 1240 
Mpc (1240 X 640 x 370 Mpc). Until recently, this was the 
largest known structur e in the Universe. U sing a friend of 
friend (FoF) algorithm lEinasto et al.l ll2014ll found that the 
linking length I = 70 h~^Mpc three syst ems from their QSR 
catal ogue coincide with the LQGs from iGlowes et ahl (l2012l . 

| 201 ^ _ 

iHorvath. Hakkila fc Bagolvl (l2013ll : 


iHorvath. Hakkila fc Bagolv ( 2014l ~) announced the dis¬ 
covery of a larger Universal structure than the Huge-LQG 
by analyzing the spatial distribution of gamma-ray bursts 
(GRBs). The 3000 Mpc size of this structure exceeds the 
size of the Sloan Great Wall by a factor of about six; This 
is currently the largest known universal structure. 

Unlike quasars, GRBs are short-lived cosmic transients 
spanning milliseconds to hundreds of seconds (see the review 
paper byll eszaro j [ 2 OO 6 I I . Due to their immense luminosi¬ 
ties, GRBs can be observed at very large cosmological dis¬ 
tances. Their host s are typically metal poor galaxies of in¬ 
termediate mass j[SayagIio ^_G4agebrook_fcj^Borgnel l2009l : 


ICastro Ceron et ah 2010l : Levesque et al. 2010l L rather than 
the massive galaxies in which quasars are generally found. 
Both quasars and GRBs should map the underlying dis¬ 
tribution of universal matter, although the details of their 
spatial distributions are not necessarily identical. Although 
the number of known GRBs for which distances have been 
accurately measured is significantly fewer than the number 
of known quasars, the surveying techniques used to identify 
these objects are more homogeneous and better-suited for 
studying structures of large angular size than are quasars. 

The existence of an object with a size of sev¬ 
eral gigaparsecs introduces questions concerning the ho¬ 
mogenous and isotropic nature of cosmological models. 
The great importance of this question necessitates fur¬ 
ther independent study into the use of GRBs for map¬ 
ping large-scale universal structures. Our analysis con¬ 
siders the space distribution of a GRB sample having 
known redshifts for the presence of large-scale anisotropies. 
The sample we use for this study is available at 
http://www.astro.caltech.edu/grbox/grbox.php As of 
October 2013, the redshifts of 361 GRBs have been deter¬ 
mined. 


2 DISTRIBUTION OF GRBS IN {R, 6, p} SPACE 

According to the cosmological principle (CP), Uni versal 
large -scale structure is homogeneous and isotropic llEllisI 
I 1975 I L The WMAP and Planck experiments have revealed 
that the Big Bang had these properties in its early expan¬ 
sion phase. The observable Universe, however, shows com¬ 


plex structures even on very large scales. The problem is to 
find a limiting scale at which the CP is valid. 

A number of well-known studies have attempted to 
fi nd the largest scale on wh ich CP is valid. According 
to lEinasto fc GramannI (ll99E ) the availabl e data suggested 
values r{t) = 130 ± 10 h~^. Yahata et al.l ll2005ll reported 
on the first result from the clustering analysis of SDSS 
quasars: the bump in the power spectrum due to the 
baryon density was not clearly visible, and they concluded 
that the galaxy distribution w as homogeneou s on s cales 
larger than 60 — lOhT^ Mpc. iTegmark et al.l (I 2 OO 6 I I us¬ 
ing luminous red galaxies (LRGs) in the Sloan Digital Sky 
Survey (SDS S) improved the evidence f or spa tial flatness 
(Dtot = 1). iLiivamagl. Tempel fc Saan (I 2 OI 2 I I have con¬ 
structed a set of supercluster catalogues for the galaxies 
from the SDSS survey main and LRG flux-limited samples. 
iBagla. Yadav fc Seshadiil 1 I 2 OO 8 II showed that in the concor¬ 
dance model, the fractal dimension makes a rapid transi- 
tion to value s close to 3 at scales between 40 and 100 Mpc. 
ISarkar et al.1 1 I 2 OO 9 II found the galaxy distribution to be ho- 
mogeneous at length-scal e s gre ater than 70 h~^ Mpc, and 
lYadav. Bagla fc Khandail (I 2 OI 0 II estimated the upper limit 
to the scale of homogen eity as being cl o se to 260Mpc 
for the ACDM model. ISdchting et al.1 (I 2 OI 2 I I studied the 
Ultra Deep Catalogue of Galaxy Structures; the cluster cat¬ 
alogue contains 1780 structures covering the redshift range 
0.2 < « < 3.0, spanning three orders of magnitude in lumi¬ 
nosity (10® < I/4 < 5 X lO^^L©) and richness from eight to 
hundreds of galaxies. These results supported the validity of 
CP. 

Assuming the validity of CP, a homogeneous isotropic 
model and a standard ACDM cosmology (Da = 0.7, Dm = 
0.3, h = 0.7, representing an Euclidean space with Da -I- 
Dm = 0.7 + 0.3 = 1) the line element in the 4D space-time 
is given by 

dP = R(t)^{dr^ + + r^sind^dip^) — dt^ (1) 

The variables in the equation have their conventional mean¬ 
ing. The change of the spatial part of dP line element in 
course of time is given by the R{t) scale factor so the spatial 
distance in the brackets i.e. 


ds^ = dr^ -I- r^dd^ + sind^dip^ 


( 2 ) 


is independent of the time. Any event in the 4D space-time 
has a footprint in the {r,6,ip\ space where r can be com¬ 
puted from the observed redshift 2 and the angular coordi¬ 
nates are given by the observations. In astronomy the an¬ 
gular coordinates are usually concretized in equatorial or 
Galactic systems. The r distance is measured by the comov¬ 
ing distance defined in the Euclidean case by 


r(z) = 


dz' 


-^0 Jq Y^r2M(l + z'^) + 


(3) 


where c is the speed of light and Ho is the Hubble constant. 

The distribution of GRBs in the {r, 0, ip} coordinate 
system can be constructed by assuming some universal for¬ 
mation history, along with spatial homogeneity and isotropy. 
This theoretical distribution, however, cannot be observed 
directly because the observations are biased by selection ef¬ 
fects. There are several factors influencing the probability 
that a GRB is detected. The limit of the instrumental sensi¬ 
tivity is such a factor, as GRBs below this threshold cannot 
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be detected. The probability of detecting a GRB depends 
also on the observational strategy of the satellite. 

The GRBs in onr sample having known redshift were 
detected by different satellites using different observational 
strategies. The method of observation results in different de¬ 
tection probabilities (known as the exposure function) since 
each satellite spends different time durations observing var¬ 
ious parts of the sky. In principle this exposure function 
can be reconstructed from the log of observations made by 
each satellite. The GRB redshift can be obtained either from 
spectral observation of the GRB afterglow or of the host 
galaxy if it can be localized. The exposure is a function of 
GRB brightness and not of GRB redshift. 

Redshift measurements have their own biases. These de¬ 
pend on the optical brightness of the afterglow or the host 
galaxy, depending on which is observed. The most significant 
factor influencing the possibility of observation is extinction 
caused by the Galactic foreground. This bias seriously in¬ 
fluences the probability of redshift measurement. However, 
galactic extinction has a known distribution and can be es¬ 
timated at cosmological distances for any angular position 
on the sky. 

Most GRB afterglows and host galaxies are optically 
faint so one needs large-aperture telescopes in order to mea¬ 
sure GRB redshifts. Northern and Southern hemisphere tele¬ 
scopes used to make these measurements are located at mid¬ 
latitudes; consequently the chance of getting the necessary 
observing time is higher in the winter than in the summer 
since the night is much longer during this season. That part 
of the sky which is accessible in the winter season has a 
higher probability for determining a GRB redshift. This bias 
is equalized when observations are made over many observ¬ 
ing seasons. However, since host galaxy measurements can 
be made at any time after the GRB is observed, this effect is 
less important when measuring the redshift from the host. 

The net effect of all these factors can be expressed by 
the following formula: 

fobs{i9, V>, z) = (p)fi„t{z) (4) 

The Ttei telescope time and Exp exposure function factors 
do not depend on the redshift. As to the Ext extinction, 
however, one may have some concerns. If the observable op¬ 
tical brightness depends on the distance, which seems to be 
a quite reasonable assumption, then the measured 2 values 
close to the Galactic equator could be systematically lower. 
Our sample of 361 GRBs, however, does not show this effect. 
One may assume, therefore, that the observed distribution 
of GRBs in the {r, 0, ip} space can be written in the form of 

fobs{'&,V’,z) = g{'d,(p)fi„t{z) (5) 

In deriving this equality it is assumed that /i„t does not 
depend on the g? angular coordinates, due to isotropy. 


3 TESTING THE ANGULAR ISOTROPY OF 
THE GRB DISTRIBUTION 

A number of approaches have been developed for studying 
the non-random departure from the homogeneous isotropic 
distribution of matter in the universe. Each method is 
sensitive to different forms of the departure of the cos¬ 
mic matter density from the homogeneous isotropic case. 


These tests have bee n previously applied to g alaxy and 
quasar distributions. iGlowes. lovino fc Shaven dlQSTl l de¬ 
scribed a three-dimensional clustering analysis of 1100 
’high-probability’ quasar ca ndidates occupying the as s igned - 
redshift band of 1.8 - 2.4. licke fc van de Weveaertl d 199 il l 
presented a geometrical model, making use the Voronoi 
foam, for the asymp totic distribution of the cosmic mass on 
10-200 Mpc scales. iGraham. Glowes fc Gampusanol dl995l l 
applied a graph theoretical method, using the minimal span¬ 
ning tree (MST), to find candidates for quasa r supe rstruc- 
tures in quasar surveys. iDoroshkevich et al.l d2004l ') used 
the MST technique to extract sets of filame nts, wall-like 
struct ures, galaxy groups, and rich clusters. IPlaten et al.l 
d201lll studied the linear Delaunay Tessellation Field Es¬ 
timator (DTFE), its higher order equivalent Natural Neigh¬ 
bour Field Estimator (NNFE), and a version of the Kriging 
interpolation. The DTFE, NNFE and Kriging approaches 
had largely similar dens it y and topology error behaviours. 

I Zhang. Springel &: Yana d201(]ll introduced a method for 
constructing galaxy density fields based on a Delaunay tes- 
sellation field estimation (DTFE). Recently. iKitaurg] ( 201,31 . 
l2014h presented the KIGEN-approach which allows for the 
first time for self-consistent phase-space reconstructions 
from galaxy redshift data. 

Even if CP is correct, the distribution of GRBs in 
the comoving frame will likely be heterogeneous due to 
the cosmic history of structure formation. However, the 
isotropy of the angular distribution can remain unchanged. 
Based on the data o btained by the CGRO BATSE exper¬ 
iment, iBriggsl (|l99.3l ~) concluded that angular distribution 
of GRBs is isotropic on a large scale. In contrast, later 
studies performed on larger samples suggest that this is 
true only for bursts of long duration (> 10 sec) rather 


than subclasses of 
sumably different 

GRBs having different durations and pre- 
progenitors llBalazs. Meszaros & HorvathI 

1998 

; iBalazs et al. 

Il99!|: iMeszaros et al.1 120001: Litvin et al.1 

2001 

;lTanvir et al. 

I 2 OO 5 I: IVavrek et al.l 2008|'). 


3.1 Testing the isotropy based on conditional 
probabilities 

The validity of Equation ([5]) can be tested by some suitably- 
chosen statistical procedure. The equation indicates that the 
distribution of the angular coordinates is independent of the 
2 redshift. By definition this independence also means that 
the g{'&,tp\z) conditional probability density of the angular 
distribution, assuming that 2 is given, does not depend on 2 , 
i.e. the joint probability density of the angular and redshift 
variables can be written in the form 

fobs{i^,(p,z) = g{i},(p\z)fir,t{z) = git9,ip)fint{z) (6) 

Let us consider an 2 i; i = 1, 2,..., n} observed sam¬ 

ple of GRB positional and redshift data. The validity of 
Equation m can be tested on this dataset in different ways. 
One approach is to test whether or not the conditional prob¬ 
ability i s independent o f the condition. This is the approach 
used bv iHorvath et akl . who split the sample into k subsam¬ 
ples according to the 2 -characteristics and tested whether or 
not the subsamples could originate from the same g{'d,ip). 
This is the null hypothesis to be tested. If the null hypoth¬ 
esis is true, then the way in which the sample is subdivided 
by 2 is not crucial. A reasonable approach is to select the 
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subsamples so that they have equal numbers of GRBs; this 
ensures the same statistical properties of the subsamples. 
The number of sub samples k is somewhat arbitrary. 

iHorvath et al.l selected k different radial bins and per¬ 
formed several tests for sample isotropy on each bin. The 
tests singled out the slice at 1.6 4 2 < 2.1 as having a sta¬ 
tistically significant sample anisotropy, which was due to a 
large cluster of GRBs in the northern Galactic sky. 

3.2 Testing the isotropy of the GRB distribution 
using joint probability factoring 

It is important to independently test the significance of the 
iHorvath et al.l result using alternate statistical methods. We 
choose to regard the GRB distribution in terms of Equation 
®. The right side of Equation m describes a factorization 
of the fobs{'&,‘P,z) joint probability density in terms of the 
angular and redshift variables. To test the validity of this 
factorization we proceed in the following way: we consider 
again the sample of Zi-,i = 1,2, ...,n} mentioned 

above. If the factorization is valid (e.g. if the angular dis¬ 
tribution is independent of redshift), then the sample is in¬ 
variant under a random resampling of the z variable while 
keeping the angular coordinates unchanged. Using this ap¬ 
proach we get a new sample that is statistically equivalent to 
the original one, assuming that the joint probability density 
factorization is valid. 

The sample coming from the fobs{'&,‘P,z) joint proba¬ 
bility density is three dimensional, so the task at hand is a 
comparison of three-dimensional samples. One way in which 
the sample can deviate from isotropy is through the pres¬ 
ence of one or more density enhancements and/or decre¬ 
ments. For this type of density perturbation fint is also de¬ 
pendent on the angular coordinates and the factorization is 
no longer valid. In this case, the resampling of the redshifts 
may change fobs as well. Comparing the original observed 
distribution with the resampled one can identify the pres¬ 
ence of density perturbations with respect to the isotropic 
case. 

To calculate the Euclidean distances in a simple way we 
introduce Xc,yc,Zc Descartes coordinates in the comoving 
r, ip frame. Using Galactic coordinates one obtains 

Xo = r cos{b)cos{l) (7) 

yo = r cos{b)sin{l) (8) 

Zc = r sin{b) (9) 

Having a sample of size N the number of objects in a dif¬ 
ferential volume dV can be written using Descartes coor¬ 
dinates in the form dN = v{x, y, z)dV = v{x, y, z)dxdydz, 
where v denotes the spatial density of the objects. It is clear 
from this formula that v{x,y,z) = Nfint{x.y.z). Obviously, 
u{x, y, z) = dN/dV. 

A trivial estimate of v is obtained by counting dN in 
a given dV. Fixing dV, the variance of u is given by dN. 
An alternative approach is to keep dN constant and look 
for the appropriate dV. This approach can be realized by 
computing the distance to the fc-th nearest neighbour in the 
sample. The distance is the conventional Euclidean one in 
our case. Proceeding in this way, u{x, y, z) = 3(A:-|-l)/(47rr^) 
where rt is the distance to the fc-th nearest neighbour. In 



k 

Figure 1. Dependence of the density variance on the fc-th order 
of the nearest neighbour. Note the rapid decrease of variance at 
lower k and the much shallower decrease for larger values. 

the following subsection we adopt this approach for the GRB 
sample. 


3.3 k-th Next neighbor statistics for the GRB 
data 

In order to find the k-th next neighbours we use the 
knn.dist(x,k) procedure in the EIVA library of the R sta¬ 
tistical packagqj. Given k, the procedure computes the dis¬ 
tances of the X sample elements up to the k-th order. Re¬ 
sampling the data r can proceed by using the sample(x,n) 
procedure of R where x refers to the dataset being resam¬ 
pled and n is the size of the new sample (which is identical 
to the original one in our case). We have made 10000 resam¬ 
plings of the comoving distance sample and computed the 
densities choosing fc = 1, 2,..., 20 for the nearest-neighbour 
distances. The densities obtained from the resampled data 
enable us to compute the mean density and its variance at 
each of the sample distances. 

At this point the value of k is quite arbitrary. Selecting a 
small value of k might make this procedure sensitive to small 
scale disturbances. However, the variance of the estimated 
density on this scale is large. In contrast, the variance of the 
density is lower for larger values of k but the density fluc¬ 
tuations of smaller scale might be smeared out. In order to 
make a reasonable compromise between k and the variance 
of the estimated density, we compute the mean variance in 
the function of k. The result is displayed in Figure [T] 

The Figure helps us to estimate the optimal selection 
of the k value. Up to fc = 5 the variance drops rapidly and 
starting about fc = 8 its change is much shallower. We select 
the values of fc = 8,10,12,14 and study densities in these 
four representative cases. The simulated 10000 runs enable 
us to calculate the mean and variance of the density at ev¬ 
ery GRB location in the sample. Using the variance and 
the mean of the density we calculate the standardized (zero 
mean and unit variance) values of the density for all points 
of the sample. In Figure [2] we display a scatter plot between 

^ http://cran.r-project.org 
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Figure 2. Dependence of the standardized logarithmic density 
on the comoving distance. The 0 (mean), la, 2a lines are marked 
with black, green and magenta colors, respectively. The GRBs 
in the Southern Galactic hemisphere and those in the Northern 
are marked with red and blue colors. Note that a group of red 
points close to the 2a line at about 2800 Mpc may correspond 
to a real density enhancement of GRBs in the Sou thern Galactic 
hemisphere. The GRB Great Wall discovered by iHorvath et al.l 
appears as a group of blue points between the Icr and 2a lines in 
the 4000-6000 Mpc distance range. 


the standardized logarithmic density and the comoving dis¬ 
tance for the cases of fc = 8,10, 12 and 14, respectively. 

A quick glance at Figure [5] re veals a group of red (South¬ 
ern) dots between the la and 2a lines at about 2800 Mpc. 
These dots may represent an as sociated group of GRBs. The 
Great GRB Wall discovered bv iHorvath et ahl can be recog¬ 
nized as an enhancement of blue (Northern) dots between 
the Icr and 2a lines in the 4000-6000 Mpc distance range. 
These impressions, however, are somewhat subjective, their 
validity needs to be confirmed by some suitable statistical 
study for a more formal significance. 

The null hypothesis in this case is the assumption that 
all spatial density enhancements of GRBs are produced 
purely by random fluctuations. Assuming the validity of the 
null hypothesis one has to compute the probability of the 
density enhancement in question. We perform a series of 
Kolmogorov-Smirnow (KS) tests and confirm that the stan¬ 
dardized logarithmic densities displayed in Figure [2] follow 
a Gaussian distribution. Denoting the logarithmic density 
by Q, the probability that the i-th element’s value is a den¬ 
sity enhancement is pi oc exp{ — g^l2) and the log-likelihood 
function has the form 

n 

L =+ const. (10) 

i=l 

The summation in Equation (HOJ yields a Xn variable with 
n degrees of freedom. All density fluctuations are restricted 
within certain distance ranges. To get a likelihood function 
sensitive to a given density enhancement and/or deficit we 


order 

00 

II 

II 

o 

k=12 

k=14 

CO 

II 

1.000 

0.914 

0.540 

0.706 

O 

rH 

II 

0.914 

1.000 

0.515 

0.661 

k=12 

0.540 

0.515 

1.000 

0.743 

k=14 

0.706 

0.661 

0.743 

1.000 


Table 1. Gorrelation between logarithmic densities computed 
from orders of nearest neighbours. 



PCI 

PC2 

PC3 

PC4 

eigen val. 

3.517 

0.222 

0.077 

0.049 

st. dev. 

1.875 

0.472 

0.279 

0.221 


Table 2. Eigenvalues resulting from the PGA. The values demon¬ 
strate clearly that only the first PG (marked in bold face) is sig¬ 
nificant. 


have to sum up successive points in Equation (IIOII . Taking 
and summing up k successive points within a distance range 
we get a x| variable having k degrees of freedom. All den¬ 
sity values are calculated from a fixed number of nearest 
neighbours in the sample. We choose for k the order of the 
nearest value from which the density is calculated. 

There is some arbitrariness to this procedure. The k 
orders of the next neighbours have been selected somewhat 
by insight. Nevertheless, the calculated densities correlate 
strongly. This property enables us to concentrate the densi¬ 
ties into one variable. 

The correlation matrix in Table [T] reveals a strong cor¬ 
relation between nearest neighbor estimation of densities of 
the orders k = 8, k = 10, k = 12 and k — 14. To get a joint 
variable we perform principal component analysis (PGA) on 
the standardized logarithmic densities used above. To get 
the PCs we use the princompO procedure from R’s stats 
library. 

By running this procedure we obtain the eigenvalues 
given in Table (2] The eigenvalues indicate the variance of 
PC’s obtained. We can infer from the variances that the 
first PC describes 91% of the total variance. We assume, 
therefore, that the information from the spatial density is 
concentrated into this variable. Prom the first PC we com¬ 
pute the Xk variables for df = 8, df = 10, df = 12 and 
df = 14 degrees of freedom. The results are displayed in 
Figure [3] 


4 DISCOVERY AND NATURE OE THE GRB 
RING 

In all frames in FigureOthere is a strong peak exceeding the 
99.9% significance level (99.95%, 99.93%, 99.96%, 99.97% at 
fc=8, 10, 12, 14, respectively) at about 2800 Mpc correspond¬ 
ing to a group of outlying points in Figure [2] This appears 
to indicate the presence of some real density enhancement. 


4.1 Discovery of the ring 

We assume that GRBs making some contribution to the 
highly significant peak shown in Figure [3] lie within the full 
width at half maximum (FWHM) angular distance of the 
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Figure 3. The calculated x\ values using the first PC. Their 
degree of freedom is given in the right top corner of the corre¬ 
sponding frame. 


GRB ID 

redshift 

distance (Mpc) 

1 (deg) 

b (deg) 

040924 

0.859 

2866 

149.05 

-42.52 

101225A 

0.847 

2836 

114.45 

-17.20 

080710 

0.845 

2831 

118.43 

-42.96 

050824 

0.828 

2786 

123.46 

-39.99 

071112C 

0.823 

2772 

150.37 

-28.43 

051022 

0.809 

2736 

106.53 

-41.28 

100816A 

0.804 

2723 

101.39 

-32.53 

120729A 

0.800 

2712 

123.85 

-12.65 

060202 

0.785 

2672 

142.92 

-20.54 


Table 3. Redshift, comoving distance, and galactic coordinates 
of the GRBs contributing to the ring-like angular structure. 


peak. Their angular distribution is shown in Figure [d] The 
most conspicuous feature in all of the frames is a ring-like 
structure in the lower left side of the frames. The redshift, 
distance and Galactic coordinates of the GRBs displaying 
the ring are given in Table |3] Using the data listed in the 
Table one can calculate the mean redshift and distance of the 
ring, along with the standard deviations of these variables, 
yielding a = 0.822, = 0.025 and = 2770 Mpc, ad = 65 

Mpc. 

By definition, the true characteristic physical size D of 
the object can be obtained from the D — Qx da = Ox 
dc/{l + z) relation, where 0 is the mean angular size and 
da is the angular distance. Substituting the corresponding 
values obtained above one gets D = 944 Mpc, corresponding 
to 1720 Mpc in the comoving frame. 



180 90 0 270 180 180 90 0 270 180 


Figure 4. Angular distribution of GRBs in the FWHM distance 
range of the highest peaks in Figure [3] The degree of freedom 
in the upper right corner has the same meaning as in Figure 0 
Note the ring-like structure of objects in the lower left part of the 
frames. 


4.2 Verification of the ring structure 

In subsection o we claim to find a regular structure in 
the shape of a ring. However, the form of this structure is 
thus far based only on a visual impression. In the following 
discussion we try to give a quantitative value supporting the 
sensibility of this subjective impression. 

The procedure we used to obtain the very low prob¬ 
ability of this density enhancement only by chance is not 
sensitive to the true shape of this clustering. Assuming this 
clustering to be real we may compute the probability of get¬ 
ting a ring-like structure only by chance by assuming some 
concrete space distribution for the objects. We make these 
calculations for the cases of a) a homogeneous sphere and 
b) a shell (for more details see the Appendix). 

The probability of observing a ring-like structure is 
p = 0.2 in the case of a shell but it is only p = 4 x 10“® 
for a homogeneous sphere. It is worth noting that real space 
distributions of cluster members generally concentrate more 
strongly towards the center than do the elements of a homo¬ 
geneous sphere. Therefore, this probability can be taken as 
an upper bound for a probability of obtaining a ring shape 
purely by chance. 

Combining this latter probability with that of observing 
the clustering purely by chance we obtain a value of p = 2 x 
10“® for observing a ring entirely by chance. Thus, despite 
the large angular size of the extended GRB cluster, we find 
evidence that the cluster represents a large extended ring in 
the 0.78 ^ 2 < 0.86 redshift range. 


4.3 The physical nature of the ring 

If we assume that the ring represents a real structure, then 
we can speculate about its nature and origins. Perhaps a 
simple explanation is that it indicates the presence of a ring¬ 
like cosmic string. This would indicate that it is a large-scale 
component of the cosmic web, representing the characteris¬ 
tic spatial distribution of the objects in the Universe. The 
main difficulty with this simple explanation lies with the 
uniformity of the redshifts (distances) along the object, in¬ 
dicating that we must be seeing the ring nearly face-on. This 
possibility cannot be excluded, but alternative explanations 
are also worth of considering. 
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GRBs are short-lived transient phenomena. The GRBs 
that compose the ring, along with their redshifts, were col¬ 
lected over a period of about ten years. The number of ob¬ 
served events is determined by the time frequency of such 
events in a given host galaxy. However, the total number of 
hosts in the region containing the GRB ring and not having 
burst events during the observation period must be much 
greater relative to those which were observed. 

The number of observed GRBs should be proportional 
to the number of progenitors in the same region, although 
the spatial stellar mass density is not necessarily propor¬ 
tional to the spatial number density of progenitors. Namely, 
the progenitors for the majority of GRBs are thought to 
be short-lived 20 — 40 Mq stars, and as such their presence 
should be strongly dependent on the star formation activity 
in their host galaxies. Thus, our knowledge of the underly¬ 
ing mass distribution is sensitive to assumptions about star 
formation within the ring galaxies. 

4.3.1 Mass of the ring 

In order to estimate the mass of the ring structure we make 
two extreme assumptions providing lower and upper bounds 
for the mass. We get a lower bound for the mass by assum¬ 
ing that the general spatial stellar mass density is the same 
in the field and in the ring’s region and only the star for¬ 
mation and consequently the GRB formation rate is higher 
here. We get an upper bound for the mass by supposing a 
strict proportionality between the stellar mass density and 
the number density of the progenitors. For both estimates 
we need to know the local stellar density. 

Several recent studies have attempted to determine the 
stellar (barionic) mass density and its re l ation to the total 
Universal mass density. [Bahcall fc Kulierl (l2014l l determined 
the stellar mass fraction and found it to be nearly constant 
on all scales above 300h~^kpc, with M*/Mtot = 0.01 ± 
0.004. iLe Fevre et all ll20ll l2014ll issued the VIMOS VLT 
Deep Survey (VVDS) final and public data release offering 
an excellent opportunity to revisit galaxy evolution. The 
VIMOS VLT Deep Survey is a comprehensive survey of the 
distant universe , cove ring all epochs since 2 ^ 6 . From this, 
iDavidzon et af.l (l2013l l measured the evolution of the galaxy 
stellar mass function from 2 = 1.3 to 2 = 0.5 using the first 
53 608 redshifts. 

iMarulli et al.l ll2013ll investigated the dependence of 
galaxy clustering on luminosity and stellar mass in the red- 
shift range 0.5 < 2 < 1.1 using the ongoing VIMOS Pub¬ 
lic Extragalactic Survey (VIPERS). B ased on their sample 
of 10095 galaxies, [Driver et al.l (l2007l l estimated the stellar 
mass densities at redshift zero amounting to 8.6 h~^ ± 0.6 x 
Mq /Mpc?. We use this local value as our measure of 
the mass density in the comoving frame and use it for our 
subsequent calculations. 

We assign a volume to the ring by computing the con¬ 
vex hull (GH) of the points representing the GRBs in the 
rest frame. Using the Qhul]0 program we obtain a value of 
1.9 X 10® Mpc® for the volume of the CH. Supposing that 
the stellar mass density is the same in the ring’s region (i.e. 

^ http://www.qhull.org/ Qhull implements the Quickhull algo¬ 
rithm for computing the convex hull 


within the GH) as in the field and only the number of pro¬ 
genitors is enhanced we compute a mass of 2.3 x 10^^ Mq 
inside the GH. 

Alternatively, if we assume that the fraction of progen¬ 
itors is the same along the ring as it is in the field, then the 
total stellar mass in the volume of a shell with 2770 Mpc 
radius and 200 Mpc thickness (the observed distance range 
of the GRBs in the ring) is 2.2 x 10^® Mq. Since the num¬ 
ber of GRBs making up the ring is about the half the total 
observed in the shell, we get a mass of 1.1 x 10^® Mq within 
the ring’s GH. 

Supposing a strict proportionality between the spatial 
densities of the number of GRB progenitors and stellar 
masses we estimate a factor of 50 times more mass than that 
which is obtained above using the field value within the GH. 
In the case of a homogeneous mass distribution within the 
CH this implies an overdensity by a factor of 50 compared 
to the field. 

In reality, however, the overdensity within the CH ap¬ 
pears to be concentrated in the outer half of its volume in 
order to produce a ring-like distribution, suggesting an over¬ 
density enhanced by a factor of more than 100. This high 
value appears to be unrealistic. To resolve this contradic¬ 
tion, the proportion of progenitors in the stellar mass has 
to be increased by at least an order of magnitude and the 
stellar mass density has to be decreased by the same factor 
in the outer half volume of the CH. This results in a value of 
1 X 10^® Mq, which still represents an overdensity of a factor 
of 10 suggesting the ring mass is in the range 10^^ —10^® Mq, 
depending on the fraction of progenitors in the stellar mass 
distribution. 

4.3.2 The case for a spheroidal structure 

To overcome the difficulty caused by the low probability of 
seeing a ring nearly face on one may assume the GRBs pop¬ 
ulate the surface of a spheroid which we see in projection. 
To demonstrate that the projection of GRBs uniformly pop¬ 
ulating the surface of the spheroid really can produce a ring 
in projection, we make MG simulations displayed in Figure 
[ 5 ] The simulations show that a ring structure can be ob¬ 
tained easily by projecting a spheroidal shell onto a plane. 
The probability of observing a ring in this way is much larger 
than that of observing a ring face on. 

Unfortunately, this approach also faces some problems. 
Assuming the observed ring is a projection onto a plane, 
one can calculate the standard deviation of distances of the 
objects to the observer. A simple calculation shows this stan¬ 
dard deviation is about 58% of the radius in case of a sphere. 
Previously we obtained 1720 Mpc for the diameter of the 
ring resulting in a 860 Mpc radius with a 499 Mpc stan¬ 
dard deviation for the comoving distances. With the pro¬ 
jection correction, however, we obtain only 65 Mpc for this 
value. This result is obviously in tension with the value of 
the standard deviation assuming a spherical distribution for 
the GRBs displaying the ring. 

The relatively low standard deviation of the distances, 
however, is not necessarily caused by some physical property 
of the structure but could be caused by the FWHM of the 
statistical signal. Increasing the distance range around the 
peak of the statistical signal relative to the value of the stan¬ 
dard deviation increases the foreground/background as well. 
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Figure 5. Monte Carlo simulation of projecting points into a 
plain, distributed uniformly on a sphere. It is worth noting that 
some of the simulations strongly resemble the observed ring. 


and indicates that the structure may be buried in the noise. 
Nevertheless, for the case of a projected sphere increasing 
the distance range in this way implies that the total number 
of GRBs displaying this structure also has to be increased by 
a factor of 6. This would cause the FWHM of the statistical 
signal to be much wider, in contrast to what is observed. 

One may resolve this tension somewhat arbitrarily in 
the following way. The 499 Mpc value for the standard devi¬ 
ation of the comoving distances was obtained by assuming a 
shell-like GRB distribution. Let us take an interval around 
the mean distance of 2770 Mpc within the standard devia¬ 
tion of 499 Mpc. The endpoints correspond to some lookback 
time difference between GRBs detected at the same moment 
by the observer. The lookback time can be calculated from 
the following equation; 

tL{z)=tH f - . (11) 

Jo (1 -t- G 

where tjy = 1/Ho is the Hubble time. Galculating the time 
difference one obtains AIl = 1.9 x 10® years. Gomputing 
the time difference taking the observed 65 Mpc standard 
deviation instead of 499 Mpc, one gets only AIl = 2.5 x 
10® years. If the GRBs displaying the observed ring really 
populate a sphere, then the presence of the low 65 Mpc 
standard deviation reveals a 2.5 x 10® year period in the 
life of the host galaxy when it is very active in producing 
GRBs. Furthermore, one has to assume that this happens 
for all hosts simultaneously. This coordinated activity may 
happen by some external effect which is responsible for the 
formation of the sphere. 


One can make a similar estimate of the sphere’s mass 
by assuming that the ring represents a real structure and is 
not simply a projection. The difference in this case is that 
the ring mass represents only a fraction of the sphere’s mass. 
There are 7 objects in the ring within the 65 Mpc standard 
deviation. There is a factor of 7 for getting the number of 
GRBs within la on the sphere. This number represents 68% 
of the total number of GRBs on the sphere, i.e. one should 
multiply the mass obtained for the ring by roughly a factor 
of ten, yielding 10^® — 10^® M0 for the sphere. 


4-3.3 Formation of the ring 


No matter whether we interpret the spatial structure of the 
ring as a torus or as a projection of a spheroidal shell, the for¬ 
mation of a structure with this large size and mass provides 
a real challenge to theoretical interpretations. In addition to 
the size and the mass of the structure, one has to explain 
why the GRB activity is much higher along the ring than it 
typically is in the field. 

There is general agreement among researchers that fol¬ 
lowing the early phase of the Big Bang the initial pertur¬ 
bations evolved into a cosmic web consisting of voids sur¬ 
rounded by string-like structures. A filamentary structure 
surrounding sp here-like voids is a typical result of gravita¬ 
tional collapse dGentrella fc MelottHl98li : Ilckelll984li . 

The hierarchy of structures in the density field inside 
voids is reflec ted by a similar hierarchy of s tructures in the 
velocity field dAragon-Galvo &: Szalavl 2013l l. The void phe¬ 
nomenon is due to the action of two processes: the synchro¬ 
nisation of density perturbations of medium and large scales, 
and t he suppression of ga laxy formation in low-density re¬ 
gions dEinasto et al.l[201lh . 

It is generally assume d that the maximal size of these 
structu res is 100—150 Mpc (|Msch_etayil99^ Emusto_^^_^ 
199^; ISuhhonenko et ^ l201ll: lAragon-Galvo fc Szalav 


2 OI 3 II . Quite recently ITuIIv et al.l d20l3 l discovered the 
local supercluster (the "Laniakea”) having a diameter 
of 320 Mpc. This scale is several times smaller than the 
estimated 1720 Mpc diameter of the GRB ring, although 
perturbations on larger s cales c annot be excluded. However, 
iDoroshkevich fc Klvpinl dl988ll have presented arguments 
that perturbations on the 200 — 300 Mpc scale should 
be excluded. This value is in a clear contradiction to 
the existence of the GRB ring and other large observed 
structures. Resolution of this contradiction is still an open 
issue. 

The existence of the ring, either as a torus or the projec¬ 
tion of a spheroidal shell, requires a higher spatial frequency 
of the progenitors along the ring than in the field. A possi¬ 
ble interpretation of the higher fraction of progenitors along 
the ring is that hosts are still in the formation process at 
6.7 X 10® years after the Big Bang. This supports the view 
that large scale structure can form and evolve slowly from 
the initial perturbations dZeldovich. Einasto fc Shandarird 
1 19821 : lEinasto et al.llioO^ I . 

Dark matter must be given a dominant role in large- 
scale structure theories in order to account quantitatively 
for the formation and evolution of the cosmic web. That 
is because the observed distributions of galaxies are incon¬ 
sistent with gravitational clustering theories and with the 
formation of super clusters in a wholly gaseous medium 
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lEinasto. Joeveer fc Saail Il980l) . Recent extensive numeri¬ 
cal studies that include dark matter indicate that its pres- 
ence accounts for the basic properties of the cosmic web 
ISpringel et al.llioosi : I Angulo et aLll2012h . and these studies 
have reproduced the cosmic star format ion history and the 
stella r mass function with some success (IVogelsberger et all 
I2OI3I 1. The very large high-resolution cosmol ogical N-body 
simu lation, the Millennium-XXL or MXXL lAngulo et all 
I2OI2I1 . which uses 303 billion particles, modeled the forma¬ 
tion of dark matter structures througho ut a 4.1 Gpc box 
in a AC DM cold dark matter cosmology. iKim et all (1201 il l 
presented two large cosmological N-body simulations, called 
Horizon Run 2 (HR2) and Horizon Run 3 (HRS), made using 
6000® = 216 billions and 7210® = 374 billion particles, span¬ 
ning a volume of (7.200 h~^Gpc)^ and (10.815 h~^Gpc)^, re¬ 
spectively. Although these sizes of the simulated volumes 
were large enough to produce very large structures, account¬ 
ing for local enhancements corresponding to the size of the 
GRB ring structure is still an open problem. We address this 
issue in the next subsection. 

4-3.4 Spatial distribution of GRBs and large scale 
structure of the Universe 

In subsection 14.3.31 we noted that very large scale cosmo¬ 
logical simulations may account for huge disturbances in 
the dark matter distri bution, in particular for LQGs an d 
the object discovered bv iHorvath. Hakkila fc Bagolvl ll2013lj . 
Nevertheless, in subsections 14.3.11 and l4l3^ we pointed out 
that the existence of the GRB ring probably can not be 
accounted for a simple enhancement of the underlying bar- 
ionic and dark matter density. Presumably, to explain the 
existence of the ring one needs a coordinated enhanced GRB 
activity in the responsible host galaxies. 

According to a widely accepted view the majority of 
the observed GRBs are resulted in collapsing high mass 
(20 — 40 Mq) stars. GRBs are very rare transient phenom¬ 
ena, consequently, they observed spatial distribution is a 
serious under sampling of the space distribution of galaxies 
in general. Furthermore, the high mass stars have short life¬ 
times, consequently GRBs prefer those galaxy hosts having 
considerable star forming activity. 

Due to their immense intrinsic brightnesses, GRBs can 
be detected at large cosmological distances. GRB 090423 
has ^ = 8.2, the lar gest spectroscopically measured redshift 
dTanvir et al.l [2009l 'l. Even though GRBs seriously under¬ 
sample the matter distribution, they are the only observed 
objects doing so for the Universe as a whole up to the dis¬ 
tance corresponding to the largest measured redshift. 

Since there is no complete observational information on 
the spatial distribution of dark and barionic matter on the 
same scale as that of GRBs, one has to use the large scale 
simulations of the distribution of the cosmic matter for mak¬ 
ing such comparisons. We used for this purpose the publicly 
available Millennium-XXL simulatior^- 

As we mentioned at the end of subsection 14.3.31 the 
4.1 Gpc size of the simulated volume is large enough to 
account for structures with characteristic size of the ring. 
Since GRBs prefer host galaxies with high star formation 

® http://galformod.mpa-garching.mpg.de/portal/mxxl.html 


Histogram of dark matter density (z=0.82) 
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Figure 6. Probability distribution of the dark matter density in 
the Millennium simulation. Distribution of dark matter density 
in the simulation (green), for galaxies in general (red), and for 
star forming galaxies (blue). The star forming galaxies prefer a 
certain range of underlying dark matter density that differs from 
that of galaxies in general. 

activity, we calibrate the dark matter density in XXL to the 
spatial number density of galaxies having large star forming 
rate (SFR) assuming that 

i's{x,y,z)) = c{Qd)Qd[x,y,z) ( 12 ) 

where Vs represents the spatial number density of star form¬ 
ing galaxies and Qd the density of the dark matter. We as¬ 
sume that the c{Qd) conversion factor depends only on Qd 
but not on the spatial coordinates and that it is identical in 
the XXL and the Millennium simulation. 

We determine the c{Qd) conversion factor using the data 
available in the Millennium simulation. The GRB ring is lo¬ 
cated in the 0.78 < 2 < 0.86 redshift range, so we select 
the 2 = 0.82 slice of the simulation. The star fo rming galax- 
ies have SFR > 30 Mq yr~^ at this redshift dPerlev et al.l 
l2015h . 

As one can infer from Figure [S] the selected star form¬ 
ing galaxies prefer a certain dark matter density range: for 
densities less than and greater this range such galaxies are 
uncommon in the sample. This range differs from that of 
galaxies in general. This may indicate that the spatial distri¬ 
bution displayed by the galaxies in general is not necessary 
identical with that shown by the GRBs. 

After determining the c{Qd) conversion factor, we obtain 
from Equation the 1/3 spatial number density of the star 
forming galaxies in the XXL simulation. Based on this spa¬ 
tial distribution we generate random samples of sizes compa¬ 
rable to that of the observed GRB frequency. For generating 
the simulated sample we use the Markov Ghain Monte Garlo 
(MGMG) method implemented in the metropO procedure 
available in R’s mcmc library. 

An important issue for using this algorithm is to check 
whether or not the simulated Markov chain has reached its 
stationary stage. We check it by computing the auto regres¬ 
sion function of the simulated sample by the acf () function 
in R. We also check the MCMG output by conventional MC. 

The known number of GRBs now exceeds a couple of 
thousand and is steadily increasing with ongoing observa¬ 
tions. Unfortunately, only a fraction of these have measured 
redshifts. Motivated by the number of known GRBs and by 
their relationship to star forming galaxies, we make MCMG 
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Figure 7. Comparison of the KS differences between the XXL 
and CSR (black) and the CSR (red) cumulative distributions of 
the k = 12*^ nearest neighbors distances. Note the significant 
differences between the XXL and CSR case at sample sizes of 
N=10000 and 20000 (bottom left and right) unlike to N=1000 and 
5000 (top left and right) where there are no significant differences. 


simulations of the Vs{x^ y, z) spatial number density of these 
galaxies from Equation getting sample sizes of 1000, 
5000, 10000 and 20000. 

We make 100 simulations for these sample sizes and 
compare them with completely spatially random (CSR) 
samples of the same sizes in the XXL volume. In subsection 
E 3 ] we computed the nearest neighbors of the fc = 8,10,12 
and 14 order. Following the same procedure here, we obtain 
the nearest neighbor distances of the k = 12^^ order for the 
XXL and the CSR samples and using the ks.testO proce¬ 
dure in R’s stats library, and compute the maximal differ¬ 
ence between the cumulative distributions. We repeat this 
procedure between the CSR nearest neighbor distributions. 
The distributions of KS differences between XXL-CSR and 
CSR-CSR samples are displayed in Figure [T] 

As one may infer from this figure, the distribution of 
KS differences between XXL-CSR samples do not differ sig¬ 
nificantly from those of CSR-CSR in the case of the sample 
sizes of Ai = 1000 and 5000. On the contrary, in the case of 
N — 10000 and 20000 the difference between the XXL-CSR 
and CSR-CSR cases is very significant. 

Based on this result one may conclude that the sim¬ 
ulated XXL samples with sizes of N — 1000 and 5000 do 
not differ from the CSR case. On the other hand, the sam¬ 
ples with sizes of N = 10000 and 20000 differ significantly 
from the CSR case. Each sample corresponds to some mean 
distance to the nearest object of the k — 12*^ order. The 
computed mean distances are tabulated in Table 01 

Obviously, groups having a characteristic size of 280 
Mpc can be detected with a sufficiently large sample size. 
This val ue corresponds to the largest structure (251 Mpc) 
found bv IPark et al.l (l2012l ~) using the HR2 simulation. The 
number of GRBs detected, however, is insufficient for reveal¬ 
ing this scale. Consequently, if the XXL simulation correctly 
represents the large scale structure of the Universe the GRBs 
reveal it as GSR on the scale corresponding to the sample 
size. 

At this point, however, it may be appropriate to repeat 


Sample size 

Mean dist. [Mpc] 

Prob. of CSR 

1000 

627 

0.39 

5000 

351 

0.71 

10000 

277 

< 2.2e- 16 

20000 

217 

< 2.2e- 16 


Table 4. Mean distances to the k = 12*^ order nearest neigh¬ 
bor of star forming galaxies at different sample sizes in the XXL 
simulation and the probability for being CSR. The significant de¬ 
viation from the CSR case is marked in bold face. Note that this 
size is consistent with the CP and more than six times smaller 
than the CRB ring in this paper. 


the remark made at the beginning of this subsection: the 
existence of the GRB ring can not be explained by a simple 
density enhancement of the underlying barionic and/or dark 
matter. It probably needs some coordinated star forming 
activity among the responsible GRB hosts. In this case the 
spatial distribution of GRBs does not necessarily trace the 
underlying matter distribution in general. However, we can 
not exclude the possibility that the XXL simulation does 
not correctly account for all possible large scale structures 
and GRBs are mapping a structure that was not simulated. 

Some concern may arise, however, concerning the in¬ 
terpretation of the ring as a true physical structure, and 
of the causal relat i onshi p between the GRBs displaying it. 
ISuhhonenko et ^ lloTJ) has pointed out that cosmic struc¬ 
tures greater than 140 Mpc in co-moving coordinates did 
not communicate with one another during the late stage 
of universal expansion preceding Recombination. The skele- 
ton of the web was created during the inflationary period 
dKofman fc ShandarinI Il988ll and evolved slowly following 
this epoch. 

The volume of the shell between 0.78 < z < 0.86 is 
20.2 Gpc^ in the co-moving frame. The corresponding vol¬ 
ume is 2.1 Gpc® for 2: = 0.2 in the case of the SDSS main 
sample and 14.4 Gpc^ for 2 = 0.4 for luminous red galax¬ 
ies (LRGs), respectively. The volume of the shell is about 
ten thousand times larger than the volume of a typ ical su¬ 
percluster found bv iLiivamagi. Temnel fc Saarl (l2012fl in the 
SDSS data. 

Since the cosmic web evolves slowly, the structure of the 
GRB ring should exhibit the same general characteristics as 
those displayed by superclusters defining the web. Compar¬ 
ing the estimated number of superclusters with the number 
of detected GRBs, it appears that every thousand superclus¬ 
ters has produced on average one measured GRB. GRBs are 
therefore very rare events superimposed on the web, and the 
small probability of GRB detection casts serious doubt on 
the nature of the GRB ring as a real physical structure. 

Taking these distributional characteristics into account 
suggests that the Ring is probably not a real physical struc¬ 
ture. Further studies will be needed to reveal whether or not 
the Ring structure could result from a low-frequency spa¬ 
tial harmonic of the large-scale matter density distribution 
and/or of universal star forming activity. 
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5 SUMMARY AND CONCLUSION 

Motivated _ by _ the rece nt discovery of 

iHorvath. Hakkila fc Bagolvl (l2013l i revealing a large 
Universal structure displayed by GRBs, we study the 
spatial distribution of these objects in the comoving frame. 
The advantage of this approach is that GRBs (which are 
short transients) have footprints in this frame that do not 
change with time. 

We assume in this approach that, for the spatially ho¬ 
mogenous and isotropic case, the joint observed distribution 
of the GRBs can be factored into two parts: one part de¬ 
pends on the angular coordinates while the other part is 
radial and depends on the redshift. 

This assumption can be tested in two different ways. 
The first method is essentially that used by iHorvath et aP 
which compares the conditional probability of the GRB an¬ 
gular distribution at different z values. The second method 
tests whether resampling the GRBs randomly makes any 
statistical changes in the distribution in the 3D comoving 
frame. 

We estimate the spatial density of GRBs by searching 
the angular separations of the k-th order nearest neighbours. 
For these computations we use the knn.dist(x,k) proce¬ 
dure in the FNN library of the R statistical package. To 
compromise between the large variance of estimated den¬ 
sities at small k values and the smearing out of real small 
scale structures at large k values, we use the spatial densities 
obtained by taking k = 8, 10, 12 and 14. 

Resampling the redshift distribution 10000 times and 
calculating the spatial densities from the samples obtained 
in this way we obtain mean densities and their variances 
assuming the null hypothesis, i.e. that the factorization of 
the spatial distribution of GRBs is valid. Subtracting the 
mean value from the observed one and dividing by the stan¬ 
dard deviation we calculate the standardized values of the 
densities obtained from the nearest neighbour procedure. 

KS tests revealed that the logarithmic densities ob¬ 
tained in this way follow a Gaussian distribution allow¬ 
ing us to get a logarithmic likelihood function as a sum of 
the squared logarithmic densities. Since the sum of squared 
Gaussian variables follows a Xk distribution with k degrees 
of freedom, by selecting objects in a certain distance separa¬ 
tion range and calculating the value of this variable we can 
test for the significance of density fluctuations. 

Since the calculated logarithmic densities in the k = 
8,10, 12,14 cases are strongly correlated pairwise, perform¬ 
ing a principal component analysis (PGA) allows us to join 
the logarithmic densities in the first (the only significant) 
PG variable representing 91% of the total variance. Com¬ 
puting Xk values from this PC for k = 8, 10,12, 14 degrees 
of freedom and plotting them as a function of the distance, 
we find a very pronounced peak at about 2800 Mpc corre¬ 
sponding to a significance of 99.95%, 99.93%, 99.96% and 
99.97%, depending on the degrees of freedom. 

We plot the angular positions of the GRBs within the 
FWHM range around the distance of the peak. Exam¬ 
ining these plots we conclude the following: 

- There is a ring consisting of 9 GRBs having a mean 
angular diameter of 36° corresponding to 1720 Mpc in the 
comoving frame. 

- The ring is located in the 0.78 < z < 0.86 redshift range 


having a standard deviation of (Jz = 0.025, corresponding to 
a comoving distance range of 2672 < dc < 2866 Mpc having 
a standard deviation of ad = 65 Mpc. 

- If one interprets the ring as a real spatial structure, 
then the observer has to see it nearly face on because of 
the small standard deviation of GRB distances around the 
object’s center. 

- The ring can be a projection of a spheroidal structure. 
Adopting this approach one has to assume that each host 
galaxy has a period of 2.5 x 10® years during which the GRB 
rate is enhanced. 

- The mass of the object responsible for the observed ring 
is estimated to be in the range of 10^^ — 10^® M0 if the true 
structure is a torus or 10^® — 10^® M0 in case of a spheroid. 

- GRBs are very rare events superimposed on the cosmic 
web identified by superclusters. Because of this, the ring is 
probably not a real physical structure. Further studies are 
needed to reveal whether or not the Ring could have been 
produced by a low-frequency spatial harmonic of the large- 
scale matter density distribution and/or of universal star 
forming activity. 
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APPENDIX A: OBSERVING A RING 
STRUCTURE BY CHANCE 

In this manuscript we have found strong evidence for a ring¬ 
like structure displayed by 9 GRBs. The probability of ob¬ 
taining this clustering only by chance is about p = 5 x 10“^, 
but this value is not sensitive to the actual pattern of the 
points within the group. Although we claim to have found 
evidence for a regular structure, the apparent shape of a 
ring is based only on a visual impression. It is useful to de¬ 
velop a quantitative measure supporting the efficacy of this 
subjective impression. 

In this appendix we compare the projection of two sim¬ 
ple spherical models; a) a homogeneous sphere and b) a 
spherical shell. It is not difficult to derive the probability 
density functions for these projections into a plain. If we de¬ 
note the projected distance from the mean position of the 
group to one of the members by q, then we can normalize 
each position relative to the maximum projection Qmax so 
that the projections vary in the range of {0,1}. In case of 
a homogeneous sphere the projected probability density is 
given by 

f{g) = - ^2 (Al) 

and in the case of a shell we get 

9 (0) = ^^ (A2) 

\/i - 

The shapes of these functions are displayed in Figure lATl 

The Figure demonstrates that the projections for a shell 
result in a significant enhancement of the points close to 
the maximal distance from the center; this is not the case 


© RAS, MNRAS 000, □-?? 





A giant ring-like structure at 0.78 < z < 0.86 displayed by GRBs 13 



Figure Al. Comparison of the projected probability densities of 
a homogeneous sphere and a shell. The median of the sphere is 
indicated. 


for a projected homogenous sphere. We can compute the 
probability of measuring all 9 points outside of the median 
distance. By definition the median splits the distribution 
into two parts of equal probability. The value of the median 
distance for a projected homogeneous sphere is Qmedian = 
0.61. 

We can calculate the probability of measuring all 9 
points in the 0.61 < q < 1 regime. To calculate the prob¬ 
ability of this case we invoke the binom.testO procedure 
available in the R statistical package. The probability of 
finding an object outside the median distance is p = 0.5 by 
definition. The probability for finding all 9 points outside 
the median is given by the binomial distribution. 

Similarly we can calculate the probability of having all 
9 objects in a region outside the Qmedian ~ 0.61 median, 
assuming that the true spatial distribution of the points is a 
shell. Integrating the g(Q) probability density in the {0.61, 
1} interval we get 


/ 9{q)= j 
0.61 0.61 


\/l “ 


= 0.7924 


(A3) 


Inserting this probability into the binomial test expression 
we get p = 0.2192. 

Summarizing the results of the tests performed above, 
we infer that the probability of observing a ring-like struc¬ 
ture from a projected 3D homogeneous density enhancement 
only by chance is p = 3.9 x 10“^ while assuming a projected 
3D shell it is much higher, p = 0.22. This gives us a good 
reason to believe that the ring is a result of a projected 3D 
shell-like regular pattern. 
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